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O ' ABSTRACT 
CO 



P-Cygni type Lya profiles exhibited in nearly half of starburst galaxies, both 
nearby and high-z, are believed to be formed by an expanding supershell 
surrounding a star-forming region. We apply the Monte Carlo code which was 



> 

o 

developed previously for static and plane-parallel medium to calculate the 
Lya line transfer in a supershell of neutral hydrogen which are expanding 

CN ! 

, radially in a spherical bulk flow. We consider typical cases that the supershell 

. has the Lya line-centre optical depth of To = 10 5 — 10 7 , a radial expansion 

velocity of V cxp ~ 300 km s 1 , and the turbulence of b ~ 40 km s 1 . We find 

o ■ 

S-^ I that there appear a few emission peaks at the frequencies corresponding to 

j§ I (2iV — l)V cxp , where the order of back scatterings N > 1. As V^ X p — * b, the 

emergent profiles become similar to those for the static medium and the peaks 
• i— i . 

■ are less prominent. We also investigate the effects of column density of the 

H ' 

' supershell on the emergent Lya profiles. We find that the number and the 

flux ratios of emission peaks are determined by interplay of r and V cxp of the 
supershell. We discuss the effects of dust extinction and the implication of our 
works in relation to recent spectroscopic observations of starburst galaxies. 

Key words: line: formation — radiative transfer — galaxies: starburst — 
galaxy: formation 



1 INTRODUCTION 

Lya emission lines render us a very useful redshift indicator for starburst galaxies at high 
redshift, because their rest-frame equivalent width was predicted to be 50 - 200 A (Chariot 
& Fall 1993). Thommes et al. (1998) suggested that the existence of strong Lya emission is a 
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necessary condition for primeval galaxies. However, dust enrichment is believed to be a fast 
process in the early evolution of galaxies, and therefore Lya should suffer severe extinction by 
dust. It has been argued that even a small amount of dust in a very thick hydrogen medium 
can destroy Lya photons because their path lengths before escape increase tremendously 
due to a huge number of resonance scatterings. Hence, dust absorption is believed to be the 
main reason for the early null detection of Lya emission from high- ,2 galaxies by narrow 
band imaging surveys. 

The Lyman break method has been very successful among a few programs to detect high- 
z galaxies. Steidel et al. (1996) presented the rest-frame UV spectra of primeval galaxies at 
(z) — 3. They continued to search for even more distant galaxies, and secured the spectra 
of galaxies at the redshift (z) = 4 (Steidel et al. 1999). The UV spectra of those primeval 
galaxies are characterized by a flat UV continuum and the three types of the Lja feature: 
symmetric emission, broad absorption in the wings, and P-Cygni type asymmetric emission. 

From HST GHRS observations of nearby H II galaxies, Kunth et al. (1998) discovered 
that P-Cygni type Lja emission lines are always accompanied by the low-ionization inter- 
stellar absorption lines blueshifted by 100 — 300 km s^ 1 with respect to the systemic velocity, 
while broad wing absorptions of Lya are accompanied by the much less blueshifted inter- 
stellar absorption lines. There are spectral similarities between nearby and high-z starburst 
galaxies. Recently there have been a number of observations to obtain rest-frame optical 
nebular lines from high-z galaxies (Teplitz et al. 2000a; Teplitz et al. 2000b; Pettini et al. 
2001). 

These observations support the idea that the kinematics of surrounding media plays a 
crucial role in the survival of Lja photons in primeval galaxies, which is similar to the case in 
nearby starburst galaxies. When the systemic redshifts were determined by using the nebular 
lines, evidence for bulk motions of several hundred km s _1 is found (Pettini et al. 2001). 
Ahn & Lee (1998) and Lee & Ahn (1998) investigated the formation of asymmetric Lya 
line profiles observed in the spectra of high-z galaxies by introducing the galactic superwind 
model. 

It is generally believed that the Lyct emission is not a good indicator for the star- 
formation rate, because of dust extinction. However, Lyct emission often bears a nearly 
unique information on high-z starburst galaxies. Therefore, we have to understand the line 
formation mechanism, in order to establish the quantitative measure of star-formation rate. 

We have been investigating the Lya line formation in thick and static media of neutral 
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Figure 1. Schematic diagram for the back-scattered photons. Here a Lya photon would propagate from Pi to P3. However, 
because the central region is devoid of neutral hydrogen, the photon traverse freely to P2, which is receding with respect to the 
photon. Details are explained in the text. 

hydrogen, and developed an accurate and efficient Monte Carlo method (Ahn, Lee, & Lee 
2000, 2001, 2002). In this paper we apply our method to an expanding spherical supershell. 
In section 2 we describe our Monte Carlo code. We show the detailed mechanism of Lya line 
transfer in expanding media in section 3, and the effects of both thickness and expansion 
velocity of the supershell in section 4. In the last section we summarize the results and 
discuss on the implication of our works. 

2 MONTE CARLO METHOD 

Our model for starburst galaxies is illustrated in Figure 1. We assume that there is a central 
star-forming region surrounded by an expanding spherical supershell. We also assume that 
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the Lyct source is the central super star cluster, because the Lya luminosity of a typical 
starburst galaxy is equivalent to ionizing photons from 10 3 — 10 4 O, B stars (Gonzalez 
Delgado 1998; Ahn 2000) and the star-forming region is compact (Frye et al. 2002). 

Expanding supershells are frequently seen in nearby starburst galaxies in a kiloparsec 
scale around the central star-forming region (Marlowe et al. 1995; Martin 1998). Recent 
HST observations revealed that the neutral column density of a supershell ranges A/hi ~ 
10 18 — 10 22 cm -2 , which is deduced from usual Voigt fitting analyses (Kunth et al. 1998). 
The Lyct line-centre optical depth is related with the H I column density Ahi by 



where b is the Doppler width for turbulence motion. Therefore, we should deal with scattering 
media with the line-centre optical depths r = 10 5 — 10 9 . However, Monte Carlo methods 
are usually inefficient in this regime, because of the large computing time due to the large 
number of resonance scatterings required before escape. Hence we developed an accelerated 
Monte Carlo method in the previous paper (Ahn, Lee, & Lee 2002). 

In an expanding supershell, some photons are scattered back into the opposite part of 
the supershell, and so they get redshifted. They experience a series of such back-scatterings 
before they escape the system. These processes enhance the escape probability of Lya reso- 
nance photons in optically thick media than in static medium, and the number of scatterings 
before escape becomes much less than that in static media. Therefore, even though we want 
to solve the problem of radiative transfer in extremely thick media, it is sufficient for us to 
use the original code described in Ahn et al. (2001). 

We illustrate the back-scattering process in Figure 1. Suppose that a photon scattered 
at Pi to a direction k || P\L% is scattered at a next scattering site P 3 if the central region 
is also uniformly filled with neutral hydrogen atoms. Since no neutral hydrogen is assumed 
to be present in P 4 P 5 , the photon traverses to P 2 , where IP5P2I = IP4P3 1 . We determine the 
distance 2d = IP4P5I along the wave vector k, and add IP4P5I to the path length IP1P3I in 
the code. Here we define the position vectors of the adjacent scatterings to be r' and r", 
r' = |r'|, and the angle 9 = ZOP1P4. An elementary geometrical calculation gives 




(1) 



d — r' cos# — Z, 



(2) 



where 



cos6> = 



r'k 



(3) 
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/ = \7\Fa\ = r'cosO - \r' 2 cos 2 9 - r' 2 + r 2 min ] 1/2 , and \Of\\ = \OP~ 5 \ = r min . 

The velocity shifts due to the bulk expansion of the medium should be taken into account. 
Referring to the geometry shown in Figure 1, the velocities v x and v 2 are given by 

vi = -Kx P cos 9 k, (4) 

and 

v 2 = Vexp ^ k, (5) 

where V cxp is the expansion velocity of the supershell and V cxp > by definition. Here s = 



\P\Pz\ is the path length corresponding to an optical depth r — — lni?, where R G [0, 1] is a 
random number. We subtract frequency corresponding to the relative velocity Av = v 2 — vi 
from the frequency x whenever a back-scattering occurs. Here we introduce the dimensionless 
frequency parameter x defined by 

--^^^ (6) 
Av D Av D y J 

which describes the frequency shift from the line-centre frequency Uq in units of the Doppler 
width Aud = vo{b/c). Here b is the Doppler parameter that combines both turbulence and 
the thermal motions of the scattering medium, and c is the speed of light. The detailed 
description of our code can be found in the previous papers (Ahn, Lee, & Lee 2000, 2001, 
2002). 

While scatterings occur in the supershell, we should also take into account the Doppler 
shift caused by the relative bulk motion between the two adjacent scattering positions. For 
the radial vectors of the adjacent scattering locations, ri and r 2 , the expansion velocities pro- 
jected into the wave vector should be considered. After a simple geometrical consideration, 
we obtain 

A " = K - (r - r) ' k ' (7) 

We subtract this shift from the frequency of a Lya photon whenever scatterings take place 
in the medium. When a Lya photon escapes the medium by a back-scattering, we consider 
the redshifts due to the final back-scattering; i.e., only vi in Eq.(4) is taken into account. We 
also calculate the path lengths of photons traversing in the supershell in order to estimate 
the amount of dust extinction. However, dust effect is considered very briefly in this paper, 
and we will concentrate on the Lya line transfer in the expanding supershell. 
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Figure 2. An emergent profile (solid line) of Ly« emission calculated by our Monte Carlo code. The dashed line represents the 
total number of scatterings before escape, and the dotted line refers to the number of successive wing scatterings just before 
escape. The horizontal axis refers to the wavelength shift in units of thermal width. The right vertical axis refers to the number 
of scatterings, and the left one refers to the number flux of the emergent spectrum. The vertical lines with numbers in the box 
indicate the frequencies of possible emission peaks. See the model parameters in the text. 

3 THE FORMATION OF P-CYGNI TYPE LYa PROFILES 

We now discuss on typical profiles of Lja resulting from the scattering processes discussed 
above. The expanding supershells in starburst galaxies have typical turbulence which is 
given by the Doppler parameter, b ~ 40 km s" 1 estimated from Voigt fittings (e.g. Kunth 
et al. 1998). Hence we assume that the Voigt parameter of the expanding supershell to 
be a = 1.5 x 10~ 4 (40 km s^/b), where the Voigt parameter a = T/4nAis D with T being 
damping constant. Here we neglect the thermal motion, because the turbulent motion usually 
exceeds the thermal motion. The FWHM of the input profiles is set to be Ax = 2, because 
the optical nebular emission lines indicate that Av = 70 ~ 100 km s _1 (Pettini et al. 2001). 
We consider a supershell with its expansion velocity being V cxp ~ 300 km s" 1 and its column 
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density 10 18 cm -2 < A/hi < 10 21 cm -2 , in accordance with observations in nearby starbursts 
(Kunth et al. 1998) and high-z galaxies (Pettini et al. 2001). These two quantities satisfy the 
condition af > 10 3 , which guarantees the validity of the analytic solution given by Adams 
(1972), Harrington (1973), and Neufeld (1990). Here we define the line centre-optical depth 
by f = y/ftTo, which is a little bit different from the previous researches. 

In Figure 2 we show the emergent profile for the model. The number of photons used 
in the Monte Carlo calculation is 80,000, and the other physical parameters are assumed 
to be the typical values mentioned above. Especially we choose the neutral column density 
A/ H i = 2x 10 20 cm~ 2 . The dotted line in the figure represents the average number of scatterings 
for a photon during its transfer, and the dashed line represents the average number of 
successive wing scatterings just before its escape. The solid line refers to the emergent Lya 
profile. The vertical bars with numbers refer to the frequency of emission peaks corresponding 
to N-th back-scattering. We can see two kinds of peaks in the profile: one located at the 
bluer frequency, and the other three at the redder part. 

We see in the figure that the weak peak at the blue part is composed of photons which 
have experienced very large number of scatterings in the wings as well as in the core. Ac- 
cording to Adams (1972), the number of scatterings for those photons is N oc r . Hence we 
can see that these are photons transferred in the supershell by so-called excursions (Adams 
1972). Some photons around x ~ are also similar ones, which is evident by the large 
number of scatterings and wing scatterings. 

These weak peaks are formed through the following processes. When photons with the 
line-centre frequency are incident upon the inner wall of the expanding supershell, its scat- 
tering is off-centred by an amount of x exp , where we define the dimensionless frequency 
corresponding to the expansion velocity x cxp = —V cxp /v t h- Hence, the fraction of transmit- 
ted photons through the supershell is approximately given by Eq.(2.23) in Neufeld (1990), 
and the profile of transmitted photons is given by Eq. (2.32) in the same paper. Since 
a^exp = —7.65 in the above simulation, we can see from Eq.(2.23) in Neufeld (1990) that a 
very small amount of photons can transmit the optically very thick supershell and majority 
of photons are back-scattered. The transmitted photons are redistributed into both sides of 
x = —x exp . The radiative transfer of those photons was studied very well both analytically 
(Neufeld 1990) and numerically (Ahn et al. 2002). 

The three peaks at the red part are composed of back-scattered photons. Since the 
supershell is receding with respect to the incident photon, scatterings usually occur at wing 
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Figure 3. Lya scattering opacity for the cases with the line-centre optical depths of tq = 10 5 , 10 6 , 10 7 , which are drawn by 
solid lines from the innermost one. The vertical dotted lines represent the effective scattering frequencies, and the number refers 
to the order of back-scattering with atoms in the supcrshcll. 

frequencies if the opacity at the frequencies is sufficiently large. Here we define the effective 
scattering opacity at a wing frequency x by T e (x). At the first scattering of Lya photon 
with the expanding supershell, the effective optical depth should be measured at x + x exp , 
i.e r e (xo + x exp )- Here we consider an input frequency xq = 0, for simplicity. Note that 
^exp < 0. If r e (xo + ^exp) ^> 1, a back-scattering occurs at the wing frequency. Then, the 
back-scattered photons are redshifted to the frequency x = x + x cxp , and they will scatter 
in the other part of the supershell along the back-scattered direction. If the supershell in 
that direction is optically thin, then those back-scattered photons can escape freely to form 
a prominent peak. In fact, that part of the supershell is also receding with respect to the 
back-scattered photons, and so the effective opacity now becomes r e (xo + 2x cxp ). The opacity 
at this frequency is even smaller than that of previous scattering, and the escape probability 
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increases. Hence some photons escape the system and form the first peak at x = Xq + x cxp . 
The photons reflected into slant directions are redshifted by an amount smaller than those 
to the normal direction, and the optical depths in those direction becomes larger. Hence 
the escape probability of those photons become smaller than that of photons in the normal 
direction. These photons contribute to the broadening of each peak. It is also noticeable 
in the figure that the primary peak is not well matched with x exp , because its blue part 
is eroded by the scattering of very thick supershell whose scattering centre is located at 
x = —x e xp in the figure. 

Lya photons repeat the above-mentioned back-scatterings until r e becomes small enough 
to permit most of the photons to escape the system. We see that photon's effective scattering 
frequency at the iV-th back-scattering is given by 2(N — l)x cxp , where N > 1. The frequency 
of photons escaping associated with iV-th back-scattering is given by (2iV — 3)x cxp . When 
N — 1, the situation is special, because the source is at rest. In this case, each frequency is 
x cxp and 0, respectively. 

We illustrate in Figure 3 how the effective opacity of successive back-scatterings varies, 
where we consider photons scattering back and forth along a diameter for simplicity. The 
effective frequencies of Lya photons successively change, as is drwan by the vertical dotted 
lines in the figure. The numbers in the box for each dotted line refers to the order of back- 
scatterings denoted by N. Let us consider the case for tq = 10 6 which corresponds to the 
case in Figure 2. We can see that the effective opacity for N — 1 is very large, and so most of 
Lya photons are back-scattered. This means that there is a little amount of transmitted flux 
with x — 0. The second or n = 2 back-scattering happens at wing frequencies much further 
away from the N — 1 frequency, and so the escape probability is enhanced, but it is still 
not large enough for most photons to escape the system. However, the emission peak forms 
at x — a; e xp, because the flux of incident Lya photons is very large. In the next turn with 
N = 3, the second peak at x — 3x cxp forms due to the process similar to the N = 2 case. 
When N = 4, we expect that the other peak, but it appears very weak because most Lya 
photons have already escaped from the system. From then on the effective opacities for the 
following back-scatterings are sufficiently small, but photons have been already exhausted. 
So no peaks at the higher frequencies appear in the emergent spectrum, as shown in Figure 2. 
We can summarize mathematically these explanation as follows. 

Ji(0) = J„ h(x exp ) ~ 0, (8) 
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h [1 - /t(^cxp)] /t(23; cxp ) 
^o/r (^ex P ) /t (2a; CX p) , 



(9) 



/3(3x ex p) 



exp ) j 



(10) 



/4(5x ex p) 



4 /r (^cx P ) /r ( 2xe xp ) / R (4x exp ) / T (6x cxp ) , 



(11) 



and so on. Generally 



ijv ([2iV - 3]x cxp ) = / /r(^cx P ) 



(12) 



N-l 



]J /r ([2n - 2]x exp ) • f T ([2N - 2]x cxp ) , 



n=2 



where the transmitted fraction (/ T ) and the reflected fraction (/ R ) are approximately de- 
scribed in Eq.(2.32) and Eq.(2.33), respectively of Neufeld (1990). Here I is the input flux, 
and I N is the flux associated with the iV-th back-scattering. 

In order to see how much emergent Lya photons are affected by dust extinction, we 
calculate the total path length of the emergent photons through the medium. Here we adopt 
iV H i ~ 2 x 10 20 cm~ 2 , b = 40 km s -1 , and V^p ~ 300 km s" 1 . We also assume that dust is 
uniformly distributed in the supershell. 

In Figure 4 we show the spectral distribution of both the total path length of Lya photons 
in the supershell (r d ) and the total number of scatterings (iV scat ). The latter represented by 
solid circles in the inset is depicted in logarithmic scale, and shows clearly two groups. One 
group is located at the upper domain with log 10 iV scat ~ 6, and they are photons having 
experienced excursions. The rest in the lower part with log 10 iV scat 1 are those having 
escaped by a number of back-scatterings. 

The crosses represent, in units of r , the total path length along which a Lya photon 
traversed in the supershell. We can see in the figure that most of emergent photons escape 
by back-scattering, which can be inferred from the small number of scatterings. It is also 
noticeable that the total path lengths are very diverse in the spectral distribution and most of 
photons traverses rather short path. Since there is no frequency-dependence in path lengths 
except for a weak trend of monotonous increases, the redder peaks in emergent profiles from 
dusty media can be destroyed at most a little bit more than the bluer peaks. This is puzzling 
because all the observed Lya profiles shows single peak which must be the N = 2 peak in 
our calculations. The solution of this puzzle may lie in the spatial distribution of dust in the 
system. If the inner radius of the dust shell is is smaller than that of neutral hydrogen, the 
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Figure 4. Spectral distribution of the number of scatterings (iV SC at) denoted by solid dots and the total path lengths (r<i) 
represented by crosses. The calculation is performed for a model with the line-centre optical depth to = 8. .3 x 10°, the expansion 
velocity in units of thermal width rr oxp = —7.65, and the FWHM of input profile Ax = 2. 

dust extinction of the higher order peak can increase. We are now undertaking this subject, 
and the results will be published in the forthcoming paper. 



4 DEPENDENCE ON THE PHYSICAL PARAMETERS 

The emergent profiles depends on a few parameters of scattering medium in our model. First 
of all, the turbulence is out of our interest, because they are constrained by observations and 
the emergent profile is dominantly formed by kinematics of the supershell in our case. Hence 
here we investigate how the Lya peaks depend on the optical depths and the expansion 
velocity of the supershell. 
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Figure 5. Emergent Lya profiles for three cases in which the supershclls have different line-centre optical depths but their 
expansion velocities are fixed. The parameters are given in the box, and the arrows refers to the location of possible emission 
peaks formed by back-scatterings. 

4.1 Effects of supershell's opacity 

We present in Figure 5 the emergent profiles for a few cases with different line-centre optical 
depths, where the expansion velocity of the supershell is fixed to be x exp = —7.65. The 
FWHM of input profiles is also set to be Ax = 2, and the Doppler parameter of the supershell 
is set to be b = 40 km s -1 . The supershell is assumed to have no dust, and the number of 
input photons are set to be equal for each calculation or N ph = 80, 000. The arrows in 
the figure represent the peaks corresponding to x exp , 3x cxp , and 5x cxp . When the opacity 
of the supershell is not sufficiently large (JV H i = 2 x 10 18 , 2 x 10 19 ), the input photons are 
transmitted without scatterings and the peak at x ~ forms, whose FWHM increases by 
a small amount compared to that of input profile. We can see in the figure that the back- 
scattered photons form the primary peak at x ~ :r ex p, the secondary peak at x ~ 3x exp , and 
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Figure 6. Emergent profiles of Lya for several cases with different expansion velocities. The horizontal axis refers to the 
wavelength in units of thermal width, and the vertical axis represents the number flux of emergent profiles. We assume that 
the input profiles have the same Gaussian input profile with Ax = 2 FWHM and the same number flux of photons or 80,000. 
The column density of neutral hydrogen in the supershell is also fixed or Nm = 2 X 10 20 cm -2 . The meaning of each line is 
shown in the box, and the arrows represent the frequencies for the primary, secondary, and tertiary peaks for x cxp = —7.65. 

the tertiary peak at x ~ 5rr cxp . As the opacity of the supershell decreases, the peaks at higher 
frequencies disappear accordingly. We have explained the mechanism of those variations in 
the previous section. 

4.2 Effects of Kinematics 

We also investigate the effects of the expansion velocity of the supershell. The results of 
Monte Carlo calculations are shown in Figure 6. We can see that the emission peaks gradually 
merge into a single peak as |x e xp| decreases. This is because the supershell behaves like a 
spherical medium which is static with respect to the Lya source as |x e xp| — *• 0. In an 
extremum with |a; C xp| < 1, the emergent profiles will have symmetric double peaks which is 
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same to those given in Ahn et al. (2002). If |rr cxp | is much larger than unit, then we can see 
the emission peaks of back-scattered photons in the emergent profile. However, as |x cxp | — > 1, 
the frequency gaps between those peaks become close, and they eventually mixed up to be 
a single broad peak. We can also expect a heavy extinction for small | 

It is also noticeable that the similar problem for the contraction supershell has a exact 
mirror-symmetry to the expanding supershell. Hence it is evident that we obtain the profiles 
to the origin x = for the contracting case. 

5 SUMMARY 

We have investigated the formation of P-Cygni type Lya by a supershell surrounding a 
star-forming region and outflowing in a bulk manner with expansion velocities of V cxp ~ 
100 — 300 km s -1 . We have modified the Monte Carlo code developed in the previous paper 
(Ahn et al. 2001), and studied in detail the Lya line transfer in the expanding supershell 
consists of neutral hydrogen gas. The Lya line-centre optical depth of the supershell is very 
large and of order r ~ 10 5 — 10 7 . Its turbulence motion is assumed to be about 40 km s _1 
in accordance with observations. 

We found that the back-scatterings play a crucial role in the formation of asymmetric 
Lya profiles. Most of Lya photons escape the system efficiently through back-scatterings. 
Whenever it is scattered back by the expanding supershell, photon's frequency is redshifted 
successively. This process promotes the redshift of Lya photons. It is also naturally expected 
that the collapsing supershell promotes blueshift of Lya photons, and the emergent profile 
for that case is simply mirror-symmetric to that of the expanding case. The emission peak 
formed by the N-th back-scattering with the supershell has a frequency of (2 A — 3)x exp , 
where N > 1. The number of peaks and their relative flux ratios in emergent profiles 
are determined by the interplay between the line-centre optical depths and the expansion 
velocities of supershells. 

Recently optical nebular lines of about a dozen of Lyman break galaxies have been 
obtained by using near-IR spectrographs and large telescopes (Pettini et al. 2001; Teplitz et 
al. 2000a; Teplitz et al. 2000b). The Lya spectra of a large number of primeval galaxies can be 
found in the literature (Steidel et al. 1996; Steidel et al. 1999; Lowenthal et al. 1997; Frye et 
al. 2002), which are reviewed in Stern & Spinrad (1999). Their UV spectra and Lya emission 
lines show similar characteristics to those of nearby starburst galaxies. Their P-Cygni type 
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Lya emission lines are always accompanied by nebular lines blueshifted with respect to 
Lya. We have showed that these asymmetric Lya can be formed by back-scatterings of Lya 
photons by an expanding supershell. 

The quantitative measurement of star-formation rates of starburst galaxies is one of the 
most crucial issues in modern cosmography. It is suggested that the expansion of surrounding 
media can help Lya photons escape the star-forming regions more freely. In light of our works, 
it is very interesting to see that nearly all the observed P-Cygni type Lya emissions from 
starburst galaxies do not show any higher order peaks which are formed by back-scatterings. 
This can certainly be attributed to dust extinction. However, it is very subtle problem if we 
want to quantify the escape fraction of Lya photons for the case of the dusty and expanding 
supershell. Hence, we believe that very careful consideration of both abundance and spatial 
distribution of dust is needed to establish a realistic model describing the star-forming 
activity in starburst galaxies whether it is nearby or primeval. We are now undertaking this 
topic, and its results will be reported in a forthcoming paper. 
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